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Abstract 

For complex potentials V(x) = — (ix) 3 — (3 2 (ix)~ 2 — 2 (3 5 (ix) 1 / 2 which are VT sym- 
metric we show that in the (3 3> 1 strong coupling regime the low lying bound states 
almost coincide with harmonic oscillators whenever the spectrum remains real (this 
means, at all 5 < 5 crit icai{f3) ~ !)• 
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1 Introduction 



Radial Schrodinger equations 



' * 1 Hr) = Eiftr), V [eff) (r) = + V(r) (1) 



characterized by the large angular momenta £ appear in molecular or nuclear physics . 
The strongly repulsive phenomenological core G/r 2 of V(r) with G ^> 1 is often 
added directly to the centrifugal term in three dimensions, £(£ + 1) = £(£ + 1) + G, 
£ = 0, 1, . . . and, in the latter context, even the effective dimensions D ^ 3 in 
£ = m + (D — 3)/2, m — 0, 1, . . . may become extremely large ||. 

Similar situations have inspired the development of a few efficient £ ^> 1 ap- 
proximation techniques applicable to many particular Hermitian models (cf. ref. || 
containing many further references). Thorough tests of their reliability are amply 
available J3J. Similar studies are missing in the non- Hermitian context, and our 
present purpose is to fill the gap. 

For the sake of definiteness, let us recollect the popular anharmonic oscillator 



£{£ + !) 



Veff(r) = + uj 2 r 2 + g r N , (2) 



which has been studied in many papers at large £. The encouraging results of these 
studies (cf. giving further references) show that the effective potential possesses 
a deep minimum at certain real coordinate r > 0. Near this value, the effective 
potential is very well approximated by the exactly solvable harmonic oscillator. Thus, 
choosing ui = and N = 4 for definiteness we have 



Veff(r) = V eff (r ) + - K//(ro)(r - r ) 2 + 0(r - r )-\ r 



2 ejjy u/v u/ u/ ' u \ 2g 

where the higher-order contributions prove small. 

The latter Hermitian example is a key guide to our present study. We shall just 
ask what happens when the potential and/or coordinates cease to be real in a way 
proposed, e.g., by Bender and Boettcher ||. The deepest physical motivation of 
the similar generalizations may probably be found in field theory where the non- 
Hermitian oscillators emerge in the most natural manner (cf. concluding remarks in 
1711 or recent considerations in HI). 
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One of the first explicit studies of a non-Hermitian model (@) has been offered 
by Caliceti et al || who revealed the amazing reality of energies in the manifestly 
complex potential V with cubic anharmonicity N = 3 and with a purely imaginary 
coupling g in one dimension (i.e., still at the real coordinates r and with the par- 
ity £ = —1,0). For a long time, this result has been treated as a mere isolated 
curiosity in the literature, in spite of the well known existence of its fully natural 
analytic continuation to the whole domain of analyticity of the potential, i.e., into 
the cut complex plane of coordinates r. The details concerning the purpose of such 
a complexification have been summarized, e.g., in ref. [ [10[| . 

A quartic, N = A parallel of the cubic model emerged a few years later in con- 
nection with the puzzling coincidence (up to a sign) between energy corrections in 
two different perturbative models An explanation has been given by Buslaev 
and Grecchi |T2|] who discovered that the connection was mediated by an auxiliary 



non-Hermitian potential (0) with N = 4, negative g and the most suitable complex 
choice of the coordinate r = x — ie where x G (— oo, oo). Its purely real and discrete 
spectrum was proved to be bounded below. 

Both the above non-Hermitian anharmonic oscillator examples share the com- 
mutativity of their Hamiltonian with the product of operations VxV = —x and 
TiT = — i interpreted as parity and time reversal, respectively, 

H = VTHVT = H t . (3) 

The relevance of such a type of anti-linear symmetry became appreciated only after 
Bender and Boettcher || re-analyzed the VT symmetric potentials (0) with complex 
couplings g ~ i N . Putting uj = for simplicity they found that the numerical and 
WKB analysis supports the hypothesis that the spectrum remains discrete, bounded 
below and purely real at any real exponent N > 2. This inspired an intensive further 
research and, within the resulting conjecture of the so called VT symmetric quantum 
mechanics Oj, the validity of eq. (Of) has been interpreted as a certain analogue or 



a weaker form of Hermiticity [pT2|- [ 14|| . 



In what follows we constrain our attention to the N = 3 anharmonicity and 
pick up a one-parametric PT symmetric generalization of eq. (fj) in section |2| and 
subsection |2.1|. In subsection 1272] we then recollect a few basic formal ingredients 
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and, in particular, explain the relation between the VT symmetry and boundary 
conditions. A way in which the smallness oil/ 1 may play a crucial role is outlined 
in the subsequent section |^. In section [| we describe our main results concerning 
the closed asymptotic representation of the energies in the two separate (viz., weak- 
and strong-coupling) regimes of the non-Hermitian cubic-plus-square-root models. 
A few numerical tests and illustrations confirm our assertions in section [5] and are 
complemented by a non-numerical discussion in section A short summary follows 
in section [7[ 



2 The problem 



2.1 Characteristic example: Cubic oscillator with spikes 

The short history of the VT symmetric quantum mechanics climaxed with the recent 
work by Dorey, Dunning and Tateo (DDT, [|TI|) who succeeded in rigorously proving 
that the generalized cubic model with u = 0, viz., 



+ 



+ 1} 



ay ir + ir" 



ip(r) = Eip(r), r = x—ie, x 6 (- 



-oo, oo 



(4) 



possesses the real and discrete spectrum whenever the angular momentum i is suffi- 
ciently large, 



> max 



■{2a -7), 



(5) 



r 

The reality of the spectrum under a suitable constraint is a characteristic feature 



of many pseudo-Hermitian models |16[]. Within this class, the PT symmetric DDT 
oscillator (f|) is distinguished by the presence of two spikes near the origin. Moreover, 



the preliminary tests performed on a simplified model in ref. [0] indicate that the 
reality of the spectrum becomes spontaneously broken. At some £ m inimai(®), only 
slightly below the bound (0), at least two energy levels have to merge and form a 
complex conjugate pair. In the other words, we may tolerate the non-Hermiticity 
of the Hamiltonian H^ DDT ^ as acceptably weak in the semi-infinite interval of the 
angular momenta £ G (£ m inimah °°) where £ m inimai grows with a for a > 5/2. 

The interference between the two spikes becomes particularly interesting in the 
strong-coupling domain of a ^> 1. In this interval, the square- root dynamics may be 
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understood as weakly non- Hermit ian (in the sense of generating the real spectrum) 
if and only if the kinematical centrifugal repulsion remains also strong, i = 0(a). 
The latter observation has attracted our attention since the non-Hermiticity may 
in general worsen the feasibility of the construction of the solutions while, as we 
already noted, many of the difficulties with the construction of the bound states may 
in principle be avoided due to the presence of the small parameter l/£. 



2.2 Complex boundary conditions 

The I— dependence in the Schrodinger equations remains the same for the real and 
complex potentials. The latter case is exemplified by our eq. (f|) and illustrated 
in Figure 1 where, for simplicity, the smooth term ir 3 is completely omitted. The 
real part of the effective potential is displayed there for the complex shift e = 0.8 
and angular momentum £ = 13/4 ps iminimaii®) at the medium a = 10. The Figure 



shows how the square-root force —a ^Jir(x) dominates at the large coordinates while 
the centrifugal spike prevails near the origin. Both these spikes would become more 
pronounced closer to the real axis and vice versa. 

With the wave functions unconstrained by boundary conditions we may construct 
two independent solutions ipip (r) of our imaginary cubic Schrodinger differential 
equation (|4]) which are analytic functions of r and/or x. In a way described in more 
detail, say, in ref. [1(| we may cut the complex plane of r = r(x) from the origin 



upwards. This means that we parametrize 

r = £expiyj, fe(0,oo), tp e (Sir /2,n/2). 
Using such a polar representation of r(x) we may assign the unique meaning to the 



square root expression Jir(x). In accordance with Figure 1 the real part of this long- 
range spike is oriented upwards, Tie^/Ir > 0. This convention makes our potential 
uniquely defined. Its complex and e— dependent effective form reads 



V[eff)[r{x)] = — ^-r a ^Jir(x) + ir 3 (x), r(x) + i e = x E (—00, 00). (6) 

Due to the analyticity of this function in the cut complex plane we may freely deform 
the integration path. The spectrum E^ DDT ^ remains unchanged for all the constant 
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shifts £ > so that one may work with the Buslaev's and Grecchi's |12[ asymptotic 



boundary conditions ip(—ie ± oo) = 0. Bender and Boettcher H emphasized that a 
further generalization of these boundary conditions may be admitted and reads 

if + e (-3tt/io, +71-/10), <p„ e (-11^/10,-7^/10). (7) 

All of these boundary conditions are mutually equivalent and form an elementary 
analytic continuation of their standard special case ip(±oo) = in the wedge-shaped 
vicinity (|7j) of both the ends of the real axis. This is illustrated in Figure 2 where the 
asymptotic wedge permitted for the cubic oscillator is marked by the symbol CO. 
Its boundary also avoids the upwards-running cut which starts at r = 0. 

In the context of the textbook quantum mechanics, the removal of the origin 
r = from our considerations has the two immediate consequences. Firstly, in the 
spirit of ref. |T2[ (where more details may be found) we may always return to the 
current Hermitian radial-equation case (in more dimensions, with r e (0,oo)) by a 
suitable limiting transition (e — > in our above notation). In this sense, one need 
not change the mathematical results but discards merely "a half of the available 
solutions as "manifestly unphysical" due to the divergence of their norm in this limit. 
Thus, for example, one simply crosses out all the quasi-even states in the solvable 
example of ref. |17| . 



Secondly, there exists an alternative physical context where the presence of the 
spikes of the form £(£ + l)/r 2 is a dynamical assumption Then, the Hermitian 



quantum system usually stays to be defined on the whole real axis, r G (— oo, oo). All 
the solutions retain their physical meaning even after the limiting transition e —>■ 
which merely represents a regularization recipe. Such a regularization is also needed 



within the so called supersymmetric quantum mechanics [19 



3 The method 

As we have seen, the non-Hermit icity is most easily introduced in a Schrodinger 
equation by a downward complex shift of the coordinate. The typical consequences 
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of this complexification may be illustrated via the simplest example (Q) with the 
vanishing g = and, say, scaled-out lu — 1. This offers us a suitable guide towards 
the l/£ approximations in their non-Hermitian generalizations. 

Although the VT symmetric g = oscillator is exactly solvable and its spec- 
trum is real, it still exhibits certain unusual features Jl7|. The energies remain 
non-equidistant and have to be numbered by the integer n = 0, 1, . . . and by the 
superscript ^\ 

£( ± ) = 4n + 2±(2£+l). (8) 

Once we abbreviate [£(£ + 1)] 1//4 = A = A(£) > and assume that this quantity is 
large, we arrive at the Schrodinger equation 



+ 



A 4 



+ r z 



ip(r) = Eijj(r). 



(9) 



We may infer that the absolute minimum of V( e j/)(r) = A A /r 2 + r 2 lies at the purely 
imaginary point r = R = — i A. In the domain of our interest with the large values 
of A ^> 1, it may be re- written in a perturbative form using the new, shifted variable 
t = r - Rq, 



2A 2 + U 2 + 



^j(t + R ) =Eij(t + R ) 



(10) 



Obviously, the textbook perturbation solution of this problem is straightforward |2( 
and shows that the contribution of the corrections is asymptotically small. We get 
the following harmonic-oscillator leading-order energy estimate for the low-lying part 
of the spectrum, 

,2 \ 



E 



n 



-2A 2 + 2(2n + l) + Cl — 



n 



0,1, 



11^ 



It is worth emphasizing that our non-Hermitian equation degenerates to its Her- 
mitian harmonic oscillator approximation which does not contain any centrifugal 
barrier. 



The numerical reliability of eq. (|TT1) is documented by the third column in Table 1. 
In the light of the Table the harmonic oscillator approximation reproduces the low- 
lying part of the toy spectrum E^ with reasonable quality even at a very small 
A = ^30 w 2.34. 
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Of course, the second, quasi-odd series E^ = 2£ + ArS +s> + 3 of energies is not 
reproduced here at all. The reason is given by the error term in eq. (P"T|). As soon as 
we moved in the domain of the extremely large i — > oo, the value of the very lowest 
quasi-even excitation energy E^ = H + 3 is already comparable with our error 
estimate. The high-lying energy E^ cannot be reproduced within the framework 
of our harmonic oscillator fit. 

4 Energies E^ DDT ^ at the large £ 

The main merit of the freedom in the choice of the shift e > in our non-Hermitian 
Schrodinger equation as well as in the related boundary conditions is that we may 
let the axis r(x) pass through a minimum of the complex interaction term. A full 
analogy with the Hermitian case is achieved in this manner. The systematic search 
for all the possible complex extremes of the effective potentials is easy and may be 
based on the elementary mathematical rule 

d r V ieff) (r)\ r=R = 0. (12) 

In the vicinity of the extreme we may approximate the unsolvable potentials by their 
reduction to the solvable harmonic oscillator wells. 

4.1 Weak-coupling domain, a <C t 

As long as the value of i is assumed to be very large in our particular DDT example, 
it makes sense to abbreviate 2£(£ + l)/3 = L 5 and replace I by the alternative large 
parameter L = L(£) ^> 1. The range of a is then limited by the condition (|5|) so 
that we may re-scale a = V 6L 5 § with < 5 ~ 1. This simplifies our Schrodinger 
equation @, 

ip(r) = L 2 E^(r), r = Ly . (13) 
The exact form of the re-scaled effective potential 



- — + L 5 W(y) 
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is L— independent. This simplifies the implicit definition of the minimum/minima of 
Vr e ff\(r) which re-scales with r = Ly to an elementary formula 



9 y W(y)\ 



y=Y 







(14) 



In writing the solutions of this equation we shall distinguish between the two separate 
intervals of a. In the weak-coupling regime with the negligible |a| <C I (i.e., vanishing 
5 ~ 0), equation ( |i"4"| ) is trivial (Y 5 = —i) and implies that the extremes of Vr e ff\(r) 
are located at the five complex points numbered by k = 0, 1, 2, 3, 4, 



Rk 



1/5 



Y, 



Y, 



_l e i(-8+4fc)7r/10 



5 = 0. 



After the overall change of scale of r = L y with L = 0(£ 2 ^ 5 ) > 1 we get our effective 
potential W(y) = V( e ff)(Ly)/L 3 re-written in the form of the Taylor series, 

3/2 



W(y) 



2^ 2 



15 



(y-Y k ) 2 + (y-Y k 



(15) 



As long as Tie Y k 4 is negative if and only if k = or k = 4 while Im Y£ is non-zero 
unless k = 2, the unique absolute minimum of the effective potential W(y) exists 
and lies at the point y = Y2 = —i. This implies that in the first two orders in our 
auxiliary small parameter \jl const < 1 we get the low- lying spectrum 



E (DDT) 



5L 3 



+ 



'15L 



(2n+l), 71=0,1, 



(16) 



We see that the choice of the integration path r(x) with e = L = i R2 enables us to 
replace our differential Schrodinger eq. (|j) by its harmonic oscillator approximation. 
The asymptotic compatibility of the related boundary conditions is illustrated by 
Figure 2. 



4.2 Strong-coupling regime, a = 0(£) 

The strength a of the long-range spike is constrained by the rule (|5]). This means 
that the VT symmetry breaks down at the couplings b~ critical ~ 1 in a way described 
in ref. [13]. The choice of 5 G (0, S cr iu ca i) remains compatible with the reality 



S 



of the spectrum and characterizes the strongly spiked regime where both our free 
parameters are comparably large, a = 0(£). 

For the small non-vanishing 5 ~ we may expect that the positions Yk of the 
five weak-coupling extremes of W (y) become only slightly shifted. In the vicinity of 
the absolute minimum with k = 2 we abbreviate r = —iLq and write 

3 



W(y)=?>/(2y I )+iy i -5jQiy 



2q< 



— 6\/6q — q 



Equation ( |14|) remains quadratic in the fifth power of the re-scaling factor in R 
-iLQ, 



d 2 

6 



Q 5 + Q 
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0. 



We abbreviate Q 5 = Z and arrive at the two eligible roots. Both of them are real 
and share the obligatory weak-coupling limit Z± = 1. When 5 grows from to 1, 
they split and move to their respective strong-coupling extremes, 



\Z± = !],_„ 



3 -, Z. = 2 - 

2' 3 



In the light of our original eq. (0) which may be re-written in the form 

1 



6' 

the larger root gives the wrong sign on the right-hand side and must be discarded for 
positive 5. The solution becomes unique and its value Z = Z- < 1 decreases with 
the growing 5. Near its absolute minimum, our effective potential W(y) depends on 
5 only via the function Z = Z(S), 



W{y) 



1 



15 



+ 5Z) + 



15 



4Z 4 / 5 



(l + Z)(y + i Z^) +... 



(17) 



As long as Z l / b (l) w 0.922, this is a fairly weak dependence. The formula leads to 
our final energy estimate 



rp(DDT) 
m 

with m — 0, 1, 
and smooth. 



H + L»+\ / 15(1 + Z) LV(2 m + !) + ■■■ 



Near 5 = the deformation of the spectrum (|I6) is continuous 
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5 Numerical tests 



Let us support the idea of applicability of our formulae by their immediate compar- 
ison with exact results generated by a suitable "brute-force" numerical method. All 
the necessary calculations will be performed with the help of the discrete variable 
representation (DVR) approach of Harris, Engerholm and Gwinn |2]J . Since its use 
is rather new in the present context, let us start from its brief description. 

In the simplest case, the DVR approach can be regarded as a variational method 
for finding bounded solutions of the Schrodinger equation. As a first step, we have 
to choose a suitable set of orthonormal functions {4>k(r)}^ =0 that should be real for 
a VT symmetric case. Then we evaluate matrix elements of the Hamiltonian: 



The approximate bounded solutions and their energies are finally obtained as eigen- 
vectors and eigenvalues of the Hamiltonian matrix H. 

Since the basis set is usually formed by special functions (e.g. orthogonal polyno- 
mials), the former integrals (p0| ) can be evaluated analytically. However, the evalua- 
tion of the latter ones (eq. is tedious. Therefore in the DVR, the Hamiltonian 
matrix is expressed and approximated in another representation. We calculate ma- 
trix elements of the coordinate operator 



which mostly can be done analytically. Then we transform the obtained matrix into 
its diagonal representation A: 




(19) 



(20) 



(21) 




(22) 



X -> A = Q- l XQ . 



(23) 



The Hamiltonian matrix is now approximated in the A-representation as 



H(A) 
T(A) 
%(A) 




T(A) + V(A) , 
Q- l TQ , 



(24) 
(25) 
(26) 
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The approximation basically arises from the fact that an eigenvector of the coordinate 
matrix X should correspond to a function that is localized around the appropriate 



eigenvalue The level of the approximation is further discussed in ref. [|22 



The DVR can also be regarded as a method that allows to construct wave func- 
tions in a set of discrete points, i.e. on a grid. For each set of basis functions we 
obtain a corresponding grid, and an appropriate approximation of the kinetic en- 
ergy operator T(A). The grid can be shifted in the complex plane or even scaled 
or rotated, provided that a corresponding scaling or rotation is also performed with 
the kinetic energy operator. The choice of basis functions affects only the spacing 
between successive points. For example, with 4>k(x) ~ sinkx, the grid is nearly 
equidistant while with the Hermite polynomials we obtain a grid that is symmetrical 
and denser in the middle. 

In our calculations we used Hermite polynomials as basis functions (see e.g. ref. 



23j for a precise definition of this DVR). The grid consisted from up to 1001 points 



from the interval x G (—17.5,17.5) with e G (0.5,5). The Hamiltonian matrix was 
diagonalized with the help of the QR routine for general complex matrices from the 
EISPACK package p4| . The results were verified to be stable with respect to the 
changes of the grid. 

As study let us first contemplate the exactly solvable eq. (0) at g — 0. 

The sample of the exact energies is given in the first column of Table 1. Their purely 
numerical DVR reproduction appears in the second column and we see that the 
method is fully reliable. 

The last two columns of the Table illustrate the efficiency and practical value of 
our non-numerical estimate (|11|). We have to emphasize that even in the domain 
of the not too large £, the difference between the exact and approximate energies 
remains small, of the order of ~ 0.5% for ground state. The error term remains 
the same also for all the excited states but this should be understood as a mere 
peculiarity of our choice of the oversimplified illustrative example. 

The genuine test of the present formulae (and of their merits as well as limitations 
of validity) only appears in Table 2. The "measure of smallness" 1/L lies there, 
roughly, in between 0.55 and 0.23. This choice makes the test quite stringent again. 
Its results are very encouraging. We witness, firstly, a quick improvement of the 
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quality of the ground state, from cca 1% at L « 2 up to the four-digit precision at 
L m 4. Secondly, the poor performance of our harmonic oscillator approximation of 
excited states at L « 2 (giving an almost 100% error already for the second excited 
state) is in a sharp contrast with the L « 4 results predicting the reasonably good 
two correct digits even for the 8th excited state. 

In Table 3 we extend the scope of our test beyond the weak-coupling regime with 
| a | <C 2£ and/or 5 = 0. The dependence of the energies on the non- vanishing values 
of the parameter a = 2^Ji(£ + 1)6 ~ 20.98 5 is shown there to agree very well with 
our asymptotic prediction even at the fairly small L w 2.36. We may note that 
the precision of our approximation appears to be almost a— independent in a broad 
range of a including also the domain of the negative values which were not discussed 
here in detail as safely protected against any possible spontaneous VT symmetry 
breaking |15|. 



6 Discussion 



At the maximal 8 ~ 1, the estimate of the low lying energies 

2/5 



E (DDT) 



I)" 



.^ L 3 + ! VL(2m+l) 
6 2 v ; 



m 



0,1, 



(27) 



differs significantly from its weak-coupling counterpart flT5D. Moreover, in the light 
of the numerical experiments of ref. |T5| we may expect the end of the applicability 
of our straightforward harmonic-oscillator approximation. At the same time, the 
underlying shift of the minimum looks insignificant. We may conclude that the 
freedom in our choice of the shift e may encounter its natural limitations near and 
beyond the value of 5 = 1. 

The latter point may comparatively easily be discussed quantitatively. Returning 
to the re-scaled form of our original Schrodinger eq. (0) we may re-write this 
complex differential equation in an equivalent form on real line, 



- — + L 3 

dx 2 



_2(x/L - ir]) 2 
Let us assume now that the shift r\ 



— 5 y6i (x/L — irf) + i (x/L — irj)' : 



(f>(x) = E(p(x) 



e/L is not too small. This enables us to 
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transform the effective potential in a series in the powers of L, 

+x 2 L (JL _ + 3 ^ + i? j 0(x) = ^(x) . (28) 

The residual term can only generate 0(1) corrections to the energies and may be 
omitted as irrelevant. As a consequence, we may once more shift the coordinate line, 
x = z + % g and eliminate the redundant linear term in z via the suitable choice of g. 
The new version of the potential reads 

(36 - 5 T] 5 + 24 r] 5 )L 2 (72 + 81 5 rj 5 + 216 r] 5 — 3 5 2 r] w + 34 5 77 10 + 12 77 10 ) L 3 

8^ ^ 2r/ 2 (36-(57/ 5 + 247f) ' 

(29) 

The derivation of the subsequent rj— dependent energies will be skipped here as 
straightforward. Their role is less important in the present context but may be 
expected to grow in the perturbative context (of course, the explicit study of the 
higher-order perturbation corrections lies already beyond the scope of our present 
paper). 

Formula fl29|) depends on a free parameter rj and offers a more flexible harmonic 
oscillator fit of the effective potential in eq. (p8|). Its 77— dependence may certainly 
prove useful in numerically oriented considerations, as it still gives the approximate 
energies up to a bounded O(R) = 0(1) error term. Non-numerically, formula (|29| ) 
may be checked as reproducing exactly our previous eq. (|15|) in the weak coupling 
limit 5 = with optimal r\ = 1. 

A small decrease of 77 = 1 — A may be contemplated as a small perturbation. The 
deep, 0(L 3 ) minimum of our //—dependent effective potential does not move in the 
first order at all. Only the shape becomes narrower pushing the 0(\/L) component 
of the energies (0) upwards by the factor 1/rj w 1 + A. This indicates a certain 
variational optimality of our previous leading-order results where the shift was in 
fact not arbitrary, rj = 1 . 
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7 Summary 



In this paper we paid attention to the complex and asymptotically cubic DDT os- 
cillator (Q). Within the so called VT symmetric quantum mechanics, this oscillator 
represents one of the most characteristic examples of a non-Hermitian (or "next to 
Hermitian") Hamiltonian with real spectrum. Although the model is not solvable in 
closed form, its appeal is enhanced by the presence of the variable coupling a and 
angular momentum £. 

We started from the observation that in the majority of the physical applica- 
tions of Schrodinger equation in D dimensions the relevant values of i are usually 
small. In this context, the DDT model itself is exceptional. In the strongly spiked 
a ^> 1 regime, it is formally consistent if and only if the angular momenta i are 
very high. In this regime, the weak non-Hermiticity (supporting the real spectrum) 
requires the presence of a strong centrifugal repulsive core. We re-interpreted such 
a descriptive property of the model as its internal, formal feature. Its mathemati- 
cal consistency is enhanced by the admissibility of the complex shifts and of a VT 
symmetric deformation of the axis of coordinates. 

In this framework, our main purpose was to find a suitable technique which would 
give the approximate low lying DDT spectrum non-numerically. This effort was 
inspired by the enormous success of the so called \ ji expansions, techniques which 
proved extremely successful within the standard, Hermitian quantum mechanics. 
Our study revealed that the transition between the Hermitian and non-Hermitian 
models is entirely smooth. We discovered, in particular, that the angular momentum 
parameter \l\ 3> 1 may serve as a guide to the introduction of the suitable harmonic 
oscillator approximation of the low-lying (in our case, DDT) spectrum. 

We may conclude that the feasibility of the harmonic oscillator approximation 
(presumably, not only in our non-Hermitian model @) is encouraging. We may 
expect that in the future the more consequent and precise solution of the similar 
complex models will prove obtainable by perturbative techniques. The leading-order 
harmonic-oscillator construction will be followed by the systematically constructed 
series of corrections in a way which would parallel the l/£ expansions for the purely 
real potentials. 
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Figure captions 

Figure 1. Real parts of the spikes in eq. flU) at e — 0.8. 

Figure 2. Optimal path r(x) and asymptotic wedges: boundary CO for the cubic 
oscillator (|]) and boundary HO for its harmonic-oscillator approximant ([15]). 

Table captions 

Table 1. Energy levels of the solvable model ([)]) with I = 5. 

Table 2. Energy levels of the generalized cubic model (^) with a = and growing I. 
Table 3. The lowest three energy levels of the generalized cubic model @ with I = 10 
and various a. 
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Table 1: Energy levels of the solvable model (H) with I = 5. 



exact 


numerical 


large-£ 


difference 


solution 


solution 


approximation 




-9 


-9.00000 


-8.954 


0.046 


-5 


-5.00000 


-4.954 


0.046 


-1 


-1.00000 


-1.954 


0.046 


3 


3.00000 


3.046 


0.046 


7 


7.00000 


7.046 


0.046 


11 


11.00000 


11.046 


0.046 


13 


13.00000 






15 


15.00000 


15.046 


0.046 


17 


17.00000 






19 


19.00000 


19.046 


0.046 
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Table 2: Energy levels of the generalized cubic model (Q) with a = and growing £. 



parameters 


numerical 


large-£ 


difference 




solution 


approximation 




£=5 


-11.52191 


-11.390 


0.132 


L=1.821 


-4.56482 


-4.000 


0.565 




1.87017 


3.390 


1.520 


£=10 


-28 76552 


-28 686 


0.079 


L=2.361 


-20 59867 


-20.271 


0.328 




-12 70640 


-11.855 


0.851 




-5 11663 


-3.439 


1.677 




2.14032 


4.976 


2.836 


£=20 


-68 72646 


-68 680 


0.046 


L=3 086 


-59 24706 


-59 058 

..J '.7 • V 7 '.7 v7 


0.189 




-49 91773 


-49.435 


0.482 




-40 74589 

_I_ V / • 1 -L ~..J v7 '7 


-39.813 


0.933 




-31 73951 

■ 7 _L ■ 1 ' 7 ' 7 '.7 -L 


-30.191 


1.549 




-22 90712 


-20 569 


2.338 




-14 25769 

-L JI ■ ^ '.7 | \J \J 


-10.947 


3.311 




-5 80054 


-1.324 


4.476 




2.45491 


8.298 


5.843 


£=50 


-211 13555 


-211.113 


0.023 


L=4.427 


1 nn Ann 

-199.68009 


-199.589 


r\ r\r\ 1 

0.091 




-188.29459 


-188.065 


0.230 




-176.98040 


-176.541 


0.440 




-165.73889 


-165.017 


0.722 




-154.57149 


-153.493 


1.079 




-143.47967 


-141.969 


1.511 




-132.46494 


-130.444 


2.020 




-121.52886 


-118.920 


2.608 
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Table 3: The lowest three energy levels of the generalized cubic model (§) with t = 10 
and various a. 



a 


numerical 


large-£ 


difference 




solution 


approximation 




20 


-58.62190 


-58.535 


0.087 




-49.83626 


-49.533 


0.303 




-41.29014 


-40.531 


0.759 


10 


-43.85223 


-43.768 


0.084 




-35.40717 


-35.083 


0.324 




-27.23003 


-26.398 


0.832 





-28.76552 


-28.686 


0.079 




-20.59867 


-20.271 


0.328 




-12.70640 


-11.855 


0.851 


-10 


-13.35529 


-13.282 


0.073 




-5.40717 


-5.089 


0.318 




2.27617 


3.104 


0.828 


-20 


2.38131 


2.447 


0.066 




10.16462 


10.463 


0.298 




17.70339 


18.478 


0.775 
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